APPLICATIONS OF NASTRAN SUBSTRUCTURING 


By R. Narayanaswami 
NASA Langley Research Center 


SUMMARY 


This paper describes the application of substructuring techniques for 
two example problems, (i) a square plate and (ii) the static analysis of a 
frame-wall interaction problem in multistory structures. Presently, multi- 
point constraint forces are not retrieved in NASTRAN. A DMAP routine for 
calculating the multipoint constraint forces is also presented herein. 
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INTRODUCTION 


The use of substructuring techniques in NASTRAN is well documented 
(references 1 and 2). However, it is felt that example problems involving 
large degree of freedom (d.o.f.) systems would bring out the advantages of 
substructuring in greater detail and will be of help to the NASTRAN user 
community. This is attempted in this paper. 

There are several cases where the analyst will be interested in evalu- 
ating the multipoint constraint forces — for example, the frame- wall 
interaction problem in multistory structures or the nuclear fuel pellet- clad- 
ding problem in nuclear engineering. These forces are not presently retrieved 
in NASTRAN. A DMAP routine, based on the Lagrange multiplier technique, is 
presented herein for the calculation of multipoint forces of constraint . 

When this DMAP routine is applied for large d.o.f. problems, the computing 
effort needed is so great as to make it impracticable. The substructuring 
feature in NASTRAN overcomes this difficulty. This paper uses the substruc- 
ture partitioning and the Lagrange multiplier technique to retrieve the 
interaction forces between the shear .wall and frame of a multistory 
structure. 


DESCRIPTION OF PROCEDURE 


The details of substructure partitioning are explained in the NASTRAN 
User’s Manual (ref. 2) and will not be described here. With reference to 
static analysis, the method is briefly outlined in the following paragraphs. 
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The complete structure is divided into a number of substructures, the 
boundaries of ■which may be specified arbitrarily; however, for convenience, it 
is preferable to make structural partitioning correspond to physical parti- 
tioning. Each substructure is first analysed separately, assuming that all 
common boundaries (joints) with the adjacent substructures are completely 
fixed, (in HASTRAN, this is called the Phase I operation. ) From this 
analysis, the displacements of all interior points in each substructure with 
the adjacent substructure boundaries fixed are evaluated. These boundaries 
are then relaxed simultaneously and the boundary displacements are determined 
from the equations of equilibrium .at the boundary joints (the Phase II HASTRAN 
operation). Each substructure can now be analysed for boundary displacements. 
Adding these to the Phase I displacements, (displacements of interior points 
in each substructure with adjacent boundaries fixed) we get the final displace- 
ments. (This is achieved in HASTRAH in Phase III operation.) 

The addition of the reduced substructure boundary loads and stiffness 
matrices to obtain the total boundary load and stiffness matrix for the 
complete structure, and the partition of the boundary displacement of the 
complete structure into the boundary displacements of the separate substruc- 
tures is achieved with the aid of partitioning vectors. The partitioning 
vector for each substructure is a vector of size n X 1 where n is the total 
degrees of freedom in the a- set. The various steps in the construction of 
the partitioning vectors are explained in ref. 2. For eases. where all the 
grid points in the total structure have been numbered distinctly, the parti- 
tioning vectors can be formed as follows: 


1. ’ Arrange the grid points in the a-set, in ascending sequence. 

List the connected degrees of freedom at these grid points 
(the components of the a-set) as scalar point internal 
indices in ascending numerical sequence starting with 1. 

This gives the size n of the partitioning vector. 

2. The partitioning vector for each substructure is obtained by 
entering real I’s in all locations where the substructure under 
consideration has connection components with any other substructure. 


The formation of the partitioning vector when one substructure has 
connection with two or more substructures and when the grid point numbering 
for the total structure shows discontinuities is illustrated in Example 
Problem Humber 2. 


EXAMPLE PROBLEMS 
Problem Humber 1 


The structural problem consists of a square plate with hinged supports 
on all boundaries. The 10 x 20 model, as shown in Fig. 1, uses one-half 
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of the structure and symmetric boundary constraints on the midline in order to 
reduce the order of the problem and the band width by one-half . Because only 
the bending modes are desired, the inplane deflections and rotations normal 
to the plane are constrained. This is the same problem as that solved in 
the NASTRAN demonstration manual (ref. j). 

The model is divided into five substructures. (This is not the best 
division of the problem; however, since the purpose herein is to demonstrate 
the use of identical substructures and the second stiffness reduction in 
Phase II, no attempt is made to choose the best subdivision. ) The a-set points 
consist of (i) points on the boundaries of the substructures, 12 thru 22, 88 
thru 98, 154 thru l64, 220 thru 230 and (ii) additional points in each sub- 
structure needed to define the dynamic response (this is largely based on the 
analyst's judgement), 55? 60, 65, 121, 126, 131, 187? 192, and 197. Note 
that Phase I runs are made only for two substructures, substructure 1 and 
substructure 2 (Sub-5 is identical to Sub-1; Sub-3 and Sub-4 are identical 
to Sub-2). There are 53 a-set points with 3 d.o.f. per grid point (Total 
d.o.f. s= 159)» Applying the boundary condition y = 0 along X = 0, 7 d.o.f. 
are eliminated; applying the condition u z = 0 X = 0 along X = 10, l4 d.o.f. 
are eliminated; this leaves 138 d.o.f. in a-set. Since all the grid points 
in the boundaries are not needed for reasonably satisfactory dynamic response 
of the structure, a second stiffness reduction is done in Phase II. The grid 
points omitted are 13 thru l6, 18 thru 21, 89 thru 92, 94 thru 97? 155 thru 158, 
l6o thru 163, 221 thru 224, and 226 thru 229 (total of 32 points each of 3 
d.o.f.). There are thus only 42 d.o.f. in the final solution of the pseudo- 
structure in Phase II. The natural frequency comparisons with and without the 
second stiffness reduction of Phase II is given in Table 1. 

Problem Number 2 

This problem deals with the analysis of a multistory structure. The 
shear wall and frame are treated as separate structures and they are discretized 
and divided into substructures as shown in Figures 2 through 5 • The shear 
wall is divided into 30 substructures (three for each story). The frame is 
divided into 10 substructures. Phase I analysis is performed for 7 of the 
30 substructures of the shear wall and 2 of the 10 substructures of the 
frame (due to the repetitive geometry, it is enough if 3 substructures of 
shear wall and 1 of frame are analysed for Phase I; however, to reduce the a-set 
points, the former approach is used). 

Substructure 1 has connection points with substructures 2 and 4; sub- 
structure 4 has connection points with substructures 1, 2, 5? and 7. The grid 
point numbering for the total structure is available, even though it is not 
continuous serially. Under these conditions, the partitioning vectors for 
all the substructures can be formed as shown below. 

(l) Arrange the grid points in the boundaries in ascending 

sequence: 49-56, 105-112, l6l-l68, 217-224, 273-280, 329-356, 

385-392, 441-448, 497-504, 1749-1758, 1805-1812, 1861-1868, 

1917-1924, 1973- 1980, 2029-2036, 2085-2092, 2i4l-2i48, 2197-2204, 
2905-2907, 2912-2914, 2919-2921, 2926-2928, 2933-2935, 2940-2942, 
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2947-29^9, 2954-2956, 2961-2963, 2968-2970, 3105-5107, 3112-3114, 
3119-5121, 3126-3128, 3155-5155, 3l40-3l42, 3147-3149, 3154-3156, 
3l6l-3l63, 3168-3170. Since each point has 2 degrees of fredom 
(u and v), the components of the a- set are listed as scalar point 
internal indices in ascending numerical sequence starting with 1 
as follows: 


A- Set 
grid point 

49 

50 

51 


Scalar point 
internal index 


3 

4 

5 

6 


1750 


145 

146 

147 

148 


2905 

2906 


289 

290 

291 

292 


3105 

3106 


349 

350 

351 

352 


3170 


4of 

408 


(2) The partitioning vector for each substructure is obtained by- 

entering real l’s in all locations where the substructure under 
consideration has connection components with any other substructure. 
The partitioning vectors for 2 sample substructures is shown on the 
following pages (size of partitioning vectors 408 x 1): 


488 



Scalar point 
Internal Index 

Substructure 1 

Substructure 

1 

1.0 

1.0 

2 

1.0 

1.0 

3 

1.0 

1.0 

4 

1.0 

1.0 

5 

1.0 

1.0 

6 

1.0 

1.0 

7 

1.0 

1.0 

8 

1.0 

1.0 

9 

1.0 

1.0 

10 

1.0 

1.0 

11 

1.0 

1.0 

12 

1.0 

1.0 

13 

1.0 

1.0 

14 

1.0 

1.0 

15 

1.0 

1.0 

16 

1.0 

1.0 

17 


1.0 

18 


1.0 

19 


1.0 

20 


1.0 

21 


1.0 

22 


1.0 

23 


1.0 

24 


1.0 

25 


1.0 

26 


1.0 

27 


1.0 

28 


1.0 

29 


1.0 

30 


1.0 

31 


1.0 

32 


1.0 


33 

34 

35 

36 

37 

38 

39 

4 0 

41 

42 

43 

44 

45 

46 

47 

48 
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Scalar point Substructure 1 Substructure 4 

Internal Index 

49 

50 

51 

52 


289 

1.0 


290 

1.0 


291 

1.0 


292 

1.0 


295 

1.0 

1.0 

294 

1.0 

1.0 

295 


1.0 

296 


1.0 

297 


1.0 

298 


1.0 

299 


1.0 

500 


1.0 


501 

502 
505 
504 


408 

The a- set points for the frame are 5, 9, 15, l6, 20, 24, 27, 51, 55s 58, 
42, 46, 49, 55, 57, 60, 64, 68, 71, 75, 79, 82, 86, 90, 95, 97, 101, 104, 108, 
112. Since there are 5 d.o.f . per grid point, (u, v, and 0 Z ) the a-set com- 
ponents total 90. The partitioning vectors for the substructures of the frame 
(size 90 X l) can be formed easily. 

Since it is of interest to know how the frame and the wall acting alone 
will resist the lateral wind load, the frame and the shear wall are analysed 
separately at first. 

This example will also be used to illustrate the use of multiple level 
substructuring. The multistory structure is to be analysed for different 
first story heights of 12 ft, 15 ft, 15 ft, and 20 ft. In order that the 
entire calculations are not to be repeated, a Phase II (initial) run is made 
where substructures 4 thru 50 of the shear wall are combined into a "super- 
substructure"; so also substructures 2 thru 10 of the frame. The Phase II 
(Final) run consists of combining the first story substructures to the super- 
substructures of shear wall and frame, respectively. The data recovery of 
substructures of interest is achieved in Phase III. For a different first- stoiy 
height of the multistory structure, the Phase I run for the substructures of 
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the first story and Phase XI (Final) runs are repeated with the necessary 
Phase III runs. 

The stiffness matrices with respect to the active degrees of freedom 
of the wall and frame, respectively, are merged by means of vector of size 
498 x 1. The interaction of the wall and frame is studied using multipoint 
constraint equations; the conditions to be satisfied being (i) u and v 
displacements at corresponding points of wall and frame are equal and (ii) 0 Z 
of frame at connection points with the wall should be equal to the fictitious 0 Z 
values of the wall obtained by dividing the difference of the vertical dis- 
placements at the two ends of the left wall at each floor level by the width. 

The value of the maximum displacement for wall alone, frame alone, and 
frame -wall interaction for the 'case of the lateral wind load is given in 
Table 2. It is to be pointed out that without substructuring, each of the 
cases investigated would have involved considerably more computing effort. 

For example, the frame- wall interaction problem has a total of 3210 degrees 
of freedom. In a direct analysis of the total structure the stiffness matrix 
of 3210 x 3210 has to be decomposed; whereas in substructuring, 6 substructures 
in each of which the size of matrix does not exceed 112 x 112, 2 of size 30 x 30, 

1 of size 5^ X 54, and 1 of size 468 x 468 are solved. 

As seen from Table 2, for solving six different problems, a total time 
of about 900 sec is only needed while using substructuring techniques whereas 
for the solution of one shear-wall problem alone, about 2450 sec is needed 
without sub structuring. The total time for solving all the cases without sub- 
structuring will be exhorbitantly high (the bulk of the time is spent on decom- 
position of the large stiffness matrix). It should be mentioned that this wide 
discrepancy in time with and without substructuring is largely due to the repeti- 
tive nature of the structure geometry of this problem and also that Phase III 
runs are performed only at the portion of interest in the structure. Nonethe- 
less, time savings are bound to result, in general, with the use of substructuring. 


EVALUATION OF MULTIPOINT CONSTRAINT FORCES 


In NASTRAN, the multipoint constraint forces are not retrieved. A 
LMAP program is written here to retrieve these forces. The theoretical basis 
for this DMAP routine lies in the use of the Lagrange multiplier technique. 

From the minimum potential energy principle, we have the functional 

TTp = 1/2 /„ d V -J % 4 f ds 

where s a is the surface upon which the tractions T are prescribed. 
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The collection of multipoint constraint equations can he written in 
the form 

0 

To account for such constraints, we invoke the method of Lagrange 
multiplers, and defining the vector j \ j- of these multipliers, we have the 
augmented functional 


7T_ 


{«’[«] W - M'M * W'H'W 


After applying the first necessary conditions. 


k. I . Cl li 

c 1° h 


we have 



This equation can now he solved for { x Note that the system of 
equations, in general, is not now positive definite, and hence the unsymmetric 
decomposition routine of NASTRAN has to he used for the solution. Prom 
the stand point of units, X’s have the unit of lh/in. or in-lh/in. 
depending on whether the particular multipoint constraint equation equates 
displacements or rotations . This discloses that from a purely physical 
standpoint, the X’s represent the average value of the distributed force or 
moment needed to satisfy the multipoint constraint equation. 


For -this problem the a-set stiffness matrix for the wall is of size 
408 x 408; that for the frame is 90 X 90 J and there are 30 multipoint con- 
straint equations. Thus the augmented matrix is the size 528 x 528. An 
unsymmetric decomposition of this matrix on CDC 6600 machine with 140K (octal) 
storage will require about 18.5 minutes. Sirice this is very expensive, an 
alternative formulation is used herein. The 4o8 x 408 a-set stiffness matrix 
of the wall is reduced to 20 x 20 (retaining only the d. o.f. at each connec- 
tion point with the frame); the 90 X 90 a-set stiffness matrix of the frame 
is reduced to 30 X 30 (retaining only the d.o.f. at each connection point 
with the wall); thus, with 30 multipoint equations, the augmented matrix 
of size 80 X 80 need only he unsymmetrically decomposed. The reduction of 
a-set stiffness and load matrices and the solution of the augmented matrix 
took only about 88 sec on CDC 6600 with 140K (octal) storage. The DMAP 
package for this frame-wall interaction problem, including the stiffness 
reductions mentioned above, is given in the Appendix. 

Even though the procedure described herein for the calculation of 
multipoint constraint forces is general and can be used for problems that 
do not involve and/or necessitate sub structuring techniques, it has to be 
emphasized that for large problems, the method can be used only with sub- 
structuring. Even then, an additional stiffness reduction would considerably 
shorten the- computing effort. This is because unsymmetric decomposition 
of large matrices will involve unacceptably high computing costs. 
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CONCLUDING REMARKS 


The NASTRAN substructuring techniques have been applied for the 
solution of a static and a dynamic problem. In both problems, sub structuring 
is found to result in considerable saving of computing effort. The multiple 
level sub structuring technique, which facilitates the efficient reanalysis 
of the structure when only a portion of the structure is modifed, has .been 
applied for a frame-wall interaction problem. In NASTRAN, the multipoint 
constraint forces are not presently retrieved. A DMAP routine for 
retrieving the multipoint constraint forces has been written and has been 
successfully used in calculating the interactive forces between the frame 
and shear wall of a multistory structure. 
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APPENDIX 


SECOND S T I FFNfc S j 5 f< £ UOC r l J.v I » P rf AS £ 2 XO R ETR J'Ev/E MPC FORGES 
LAGRANGE MULTIPLIER TeZMNIJUE 


MAST K A N SOURCE P k 0 G K A M COMPILATION 
CMAP-DMAP INSTRUCTION 
NO. 


BEGIN i 
* 

* RG&A,PGA ARc ST I Ft- AND La AO MATRICES 0? WALL CSU£ 4jbX40d) 

S K.G0b,P&6 A Ac STIFF AKlO LOAD MATRICES Uf FRAME (SUE *0X90) 

$ SlNtc THEkE IS NO LUAU ON THE WALL i T HE PGA VECTOR tt-lLL NOT BE USED 
i KGi,KG2 ARE REDUCED STIFF MATRICES OF w'ALLUOX.iO ) AN0 FRAME 1 30X30) 

J, KGtf ANU PG ARE STIFF ANO LOAO MATRICES FCTK FRAME-WALL COMBINED 
i C Ma f LS Th£ MATRIX OF MULTI-POINT CONST KA 1 All EQUATIONS 

* STIF ANO LOAO ARE THc AUGMENTED STIFFNESS AND LOAD MATRICES 

$ UELTA IS VECTOR Uf u LSPLACEMfcNl S» LAMBDA IS VECTOR OF MFC FORCES 
4 

INP 0 TT 1 / , , , » / C , N » - 1 /C , N , D/C tNiFRArtVfALL 4 
INPUTTi / , , , ,/ C » N, - U G , N* 1 /C » N»USE RTP i 4 
I N Pu T T 1 /KGGA , PC A , , r/C ,N> 1 0/C, N. I 4 
INPUTTI . /KGGBvPCe, * »/C» N » 29 4 


SM Pi USeTl, KGGA »# ,/OOA, KG1 , KUDO A, lGOA » UOuA » imi 4 

5MP i USETUAGGS, , « /CUB > RG2. »KC0G6 » L OLB # UOOB » » » .» 4 

8 5S02 U5ET2 , * » »GD 8, ,P06/ , PUS, , PG2 4 

9 MERGE KOI, , » KG2t VEC T 3 , /xGG/C ,N L/C »Nt 2 /C ,N,6 ■> 

10 MERGE, , PG2 , , , , V6CT 3 / PC-/ C , N , 1/G* N , 2/C, N ,2 4 

U TRNSP CM AT/C TRNSP $ 

U MERGE kG&,CTRNSP,CMAT , , v ECT4»/ST IF/ C,N« ”1/C» N, 2/ C» N» 6 4 

13 MERGE P G , , f » , V EC T 4/ L OAU/ C , N , 1 /C , N , 2 VC , N ,2 4 

l A SOLVE ST IF , L 040/ SOLVECT/ C , N , 0/ C , N, 1 /C, N , 2 / <?, N , 2 4 

15 PARTN SOL VfcC T , » VECT4/DEL TA, LAMBDA , » /C, N , 1 /C , N, 2/C , N , i/C *N 1 2 4 

16 MATPRN Ofc“LTA, LAMSOA, >♦// 4 

17 END 4 

**N0 EKKOhS FOUNO - EXECUTE NASTRAN PROGRAM** 


SSkTM 8 U L K D 4 T A ECHO 


1 

.. z 

.. 3 

... 4 

• * 5 

.. 6 

.. 7 

a 

9 

. . 10 

DTI 

UShTl 

0 

0 

4C?S 

0 

500 

0 

0 


Oil 

UScT i 

i 

406 

4 

116 

litf 

116 

life 

+ 101 

4 L 01 

Hi 

ns 

US 

116 

1 Lo 

litf 

Utf 

litf 

+ 102 

+ 102 

iii 

U6 

496 

496 

116 

116 

utf 

litf 

+103 

+ 103 

116 

ii6 

il 6 

US 

A 16 

116 

UG 

Utf 

+ 104 

+ i £4 

li.4 

11 6 

4-56 

4*6 

i 16 

Utf 

utf 

iltf 

+ 105 

+ 10 5 

ii6 

1 16 

116 

116 

•Utf 

116 

Utf 

Utf 

+ 106 

+106 

Li 6 

116 

4 96 

496 

lit 

116 

Lit 

Utf 

♦ 10 7 

+ iGT 

116 

ii6 

116 

116 

116 

Utf. 

life 

116 

+156 

+ 10S 

1 i 6 

1 1 6 

496 

496 

life 

116 

116 

Utf 

+ 109 

+ 10 9 

116 

ii6 

116 

i!6 

Utf 

116 

Utf 

US 

+110 

♦ ilO 

1 1 6 

LlS 

4 96 

496 

u6 

lib 

Utf 

Utf 

+ 1 LI 

mi 

i i.6 

ilG 

U6 

US 

life 

116 

i 16 

Utf 

+ 112 

+ iil 

iU 

li^ 

496 

496 

Utf 

lib 

Utf 

litf 

+ 1 1 3 - 

+ 113 

116 

116 

Utf 

116 

Lib 

116 

Utf 

116 

+11 4 

+ 1 L4 

ilS 

U £ 

4 96 

496 

116 

Utf 

116 

litf 

+115 

+ 1 i 5 

iii 

Utf 

116 

lit 

lib 

116 

Utf 

1 16 

+ 116 

+ 1 16 

U6 

UG 

496 

496 

116 

116 

Utf 

116 

+ 117 

HIT 

ii6 

U6 

116 

US 

116 

Utf 

utf 

116 

♦ 116 

♦ lie 

11 5 

L 16 

li6 

116 

116 

116 

116 

116 

♦ 119 

+ U9 

lifo 

11 6 

Lit 

116 

116 

lib 

Utf 

litf 

+120 

+ 12 0 

116 

116 

U 6 

116 

life 

116 

litf 

UG 

+ 121 

+m 

116 

116 

116 

116 

no 

utf 

Utf 

1 itf 

+122 

+ 122 

i 16 

lU 

1 L6 

116 

lit 

116 

Utf 

litf 

+ 123 

♦ 125 

116 

116 

116 

116 

116 

lltf 

Utf 

iltf 

♦ 124 

+ 1'2 4 

116 

116 

u6 

116 

116 

life 

Utf 

Utf 

+125 

+12 5 

LLi 

tl6 

116 

Utf 

116 

litf 

i 16 

1 ife 

+ 12 6 

+ 126 

iiS 

lx6 

116 

116 

116 

116 

116 

116 

+ 127 

♦ 127 

li 6 

Lit 

116 

116 

Utf 

Utf 

116 

116 

+12 8 

+ Ug 

11^ 

lit 

116 

Ufl 

1 16 

lie; 

lib 

utf 

+ 129 

+ 129' 

lii 

LL6 

1 LG 

116 

116 

utf 

life 

Utf 

+ 130 

+ i3<? 

116 

LL6 

116 

U6 

1 16 

utf 

Utf 

1 Lu 

+ 131 

+ I3i 

liG 

116 

116 

116 

iio 

utf 

Utf 

116 

+13*2 

+ 132 

UG 

US 

116 

116 

Utf 

litf 

Utf 

life 

♦ 133 

*133 

iU 

116 

Lit 

UG 

life 

litf 

life 

116 

+ 134 

* 1 34 

UG 

Alt 

116 

life 

116 

litf 

116 

litf 

+13 5 

+ 135 

U6 

1 i 6 

1 16 

lit 

116 

litf 

litf 

Utf 

+ 136 

+ i 36 

1 16 

no 

116 

X lS» 

116 

116 

1 16 

UG 

♦ 13 7 

+ 13? 

116 

116 

1 16 

Aitf 

116 

Utf 

litf 

11& 

+13 8 

+ 1 Jfe 

U6 

116 

116 

Utf 

utf 

Utf 

116 

life 

+ 139 
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u€ 

116 

i 16 

u6 

Utf 

Utf 

Utf 

Utf 

+140 

+ 190 

L 16 

116 

1L6 

ii6 

116 

116 

i 16 

Utf 

+ 141 

+ 191 

116 

lib 

1 16 

utf 

116 

116 

116 

116 

+ 142 

+ 192 

lit 

il6 

116 

ns 

116 

116 

116 

116 

+ 143 

+ 193 

116 

U6 

116 

Lib 

utf 

116 

116 

116 

+ 144 

+ i4M 

Hi 

lie 

116 

UG 

il6 

life 

utf 

Utf 

+ 145 

♦ 145 

M6 

U6 

116 

U6 

116 

116 

116 

116 

+ 1 4 6 

+ 146 

Ii5 

iio 

LX6 

UG 

utf 

116 

116 

116 

+147 

+ 1+7 

116 

116 

no 

UG 

U6 

utf 

116 

116 

+ 148 

♦ 148 

116 

U6 

no 

116 

life 

Utf 

life 

life 

+ 149 

+ i4 9 

iU 

ii6 

Lib 

life 

A 16 

116 

iib 

life 

+ 150 

♦ 1 50 

116 

116 

LIS 

ii6 

i 16 

116 

life 

116 

+ 151 

+ 151 

116 

lie 
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d r i 

U5ET2 
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0 

90 

0 

500 

0 

0 ■ 


LT l 

1/StTi 

1 

U6 

116 

Lit 

lib 

LL& 

life 

♦201 

+201 

496 

496 

496 

Utf 

116 

116 

116 

iltf 

, +202 

+ 202 

116 

4 96 

4 96 

456 

1 16 

116 

life 

life 

+203 

+2p3 

116 

116 

496 

496 

496 

lit 

life 

Utf 

+204 

+ Z0 4 

ll6 

116 

116 

^56 

496 

496 

life 

iltf 

+205 

+ 205 

116 

116 

US 

116 

4 96 

496 

*»96 

116 

+206 

+206 

116 

11 6 

116 

11 G 

116 

4 96 

496 

<96 
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+20 7 

Lli 

Utf 

116 

lU 

116 

116 

4 96 

496 

+208 

+20 G 

456 

AiO 

LL6 
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Table 1. Natural Frequency Comparisons, CPS 


Mode 

Theoretical 

NASTBAN 
10 x 20 

3 Substructures 
(No stiffhess 
reduction in Fh. 2 ) 

5 Substructures 
(Stiffhess 
reduction in Fh. 2 ) 

1 

0.9069 


0.9056 


0.9063 

0.9068 

2 

2.2672 


2.2634 


2.2707 

2.2812 

3 

4.5345 


4.5329 


4.6070 

4.6629 

Total (Phase 1 and 
Phase 2 ) time on CDC 
6600 machine 

258.5 sec 

238.5 sec 


Table 2. Value of maximum displacement, in. ft., 
due to lateral wind load of 1 kip/sq. ft. 



Frame 

Alone 

Shear- vail 
Alone 

Frame-Wall interaction 
problem 

Case 1-First 
story height = 12 ' 

0.051145 

0.019336 

o«oi890i 

Case 2-First 
story height * 13 * 

0 .C 52557 

0.019855 

0.019448 

(Total time for 
Fh. 1 and Fh. 2 
on CDC 6600) 

50 + 30 
* 80 sec 

290 + 220 
* 510 sec 

160 + 140 
= 300 sec 


Total time on CDC 6600 for solving the problem of shear vail alone ) 

(with a First story height of 15 ') without the use of substructuring * 2450 seC. 

but using omit d.o.f» ' 


495 






























